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Abstract 

We construct integrable generalizations of the elliptic Calogero-Sutherland- 
Moser model of particles with spin, involving noncommutative spin interac¬ 
tions. The spin coupling potential is a modular function and, generically, 
breaks the global spin symmetry of the model down to a product of U( 1) 
phase symmetries. Previously known models are recovered as special cases. 
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1 Introduction 


The inverse-square interacting particle system ii 1 and its spin generalizations 
0, I, 0 , 0 i I are important models of many-body systems, due to their exact 
solvability and intimate connection to spin chain systems [10, [11], 12, 13, 14 , 2- 
dimensional Yang-Mills theories [[15], [16], [I7j| etc. Reference |18| is a classic review, 


while |L9], pOL pl|] present newer different perspectives. 


The prototype of these models is the spin-Calogero (‘rational) scattering model 
of particles on the line carrying U( n) spin and interacting with two-body inverse- 
square potentials with a U (n)-invariant spin coupling. Most other models can be 
obtained as appropriate reductions of this model, taking advantage of its discrete 

m 


or continuous symmetries 


In particular, generalizations involving U( n) 


non¬ 


invariant interactions can be obtained this way, recovering the trigonometric models 
derived in |23|, [24f and extending them to the elliptic case [|221. 


An unrelated development has been the recent progress in nonconnnutative field 
theory and matrix models. Spatial noncommutativity can be traced back to Heisen¬ 
berg and naturally arises in lowest Landau level physics [25]. Its current manifesta¬ 
tion originates in matrix, string and membrane theory |26|, |27| and came into focus 
with the work of Connes, Douglas and Schwartz [p| . 


So far these two fields remained unrelated. In this letter we show how they can 
be cross-fertilized by borrowing notions of nonconnnutative geometry and applying 
them in the reduction scheme of the Calogero model to obtain a new integrable 
elliptic model involving non-tZ(n) invariant nonconnnutative spin interactions. Such 
a modification of the spin interaction may serve to test the ‘flavor stiffness’ of the 
original spin model, to stress the degeneracy structure of the energy spectrum and 
to identify universality features of this class of models. 


2 The reduction scheme 


The basic technique that we will use consists of reducing a system of infinitely many 
particles with spin to a finite system with generalized interactions. The reader 
should refer to p2|, and especially to the elliptic case with spin, for a more detailed 


description of the method. 

The starting point is the spin-Calogero system with classical U( n) degrees of 
freedom. This system can be obtained, for instance, from the model in SB (which 
can itself be obtained as a reduction of a hermitian matrix model into nontrivial 
angular momentum sectors) by redistributing the global U(N ) degrees of freedom of 
this model into individual particle spins, or, alternatively, from the infinite-volume 
classical limit of the spin model derived and solved in [TIJ . The hamiltonian for N 
particles reads 


N 




tr (SiSj) 



















































































Xi and Pi are one-dimensional canonical coordinates and momenta; S t are a set of 
independent classical U (n) spins of rank one and length £, that is, n x n rank-one 
hermitian matrices satisfying 

tr(Si) 2 = f (2) 

and with Poisson brackets 


{('S'*)a6, (S'jOcd} = [{Si)ad$cb - 5ad (-%)<*] 


( 3 ) 


Such spins can be realized in terms of oscillators [|16] : 

(Si)ab = a a i A\ , a, b = 0 ... n - 1 


( 4 ) 


where (Af , A?) are a set of niV independent classical harmonic oscillator canonical 
pairs with Poisson brackets: 


{A a i ,A b j } = i8 ij 8 t 


ab 


( 5 ) 


and satisfying the constraint 


= £ for all i (6) 

a 

In the above model we can analytically continue the coordinates Xj and momenta 
Pi = Xi to the complex plane. The integrability and solvability of this model trivially 
extends to the complex case. Such extensions will be useful, provided that we can 
identify a real subsystem, which will be the physical system of interest. 

The hamiltonian H is invariant under particle permutations, coordinate space 
translations and global spin rotations: 

i^P N (i) i Xi^Xi + c] Si - USiU - 1 ( 7 ) 


where p N is an element of the permutation group SV, c is a constant complex param¬ 
eter and U is a constant U( n) matrix. We can arrange, therefore, for special initial 
configurations that are invariant under some combination of the above symmetries 
and be sure that these conditions will be preserved in time. 

We shall choose the configuration to be a replication of N real coordinates over 
an infinite complex lattice with periods ci and C 2 . The N particles effectively live 
on the real coordinate of a complex torus. The total number of particles on the 
covering space is, thus, infinite and we can parametrize them with the triplet of 
indices (i;m,n), where i = 1,.. .N labels the particles in each cell and m, n G Z 
parametrize the cell. Shifts in m and n are elements of the (infinite) permutation 
symmetry of the system. The kinematical variables are chosen to obey 


i;m+l,n T Cj , Pi\m-\-l,n Pv,m,n 

1 T C 2 , Pi\m,n +1 Pi-,m,n 

ensuring that we are dealing with lattice copies. This means 

T 771C\ XIC 2 , Pi\m,n Pi 


( 8 ) 

( 9 ) 

( 10 ) 
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The above conditions, being an invariance under combined permutations and trans¬ 
lations, are dynamically preserved. To also preserve the condition that X{ are real, 
we should impose invariance under the imaginary parity transformation x —> x*. 
Choosing ci real, the only possibilities for c 2 are 


C 2 + c 2 = 0 or c 2 + c 2 = Ci 


( 11 ) 


The first choice (c 2 imaginary) leads to an orthogonal lattice, while the second choice 
leads to a rhombic lattice. 

To ensure full preservation of the lattice structure, we should also impose appro¬ 
priate periodic conditions for the spins. In this, we can take advantage of the global 
spin symmetry of the system and impose 


*S);m+l ,n 

*Si;m,n+l 


TTQ. 

T/<v 
v u i ;m,n 


u 

V 


-1 

-1 


( 12 ) 


with U and V two constant matrices. That is, spins can pick up U( n) transforma¬ 
tions as they move around the cycles ci or c 2 . Consistency requires that <S'i ;m +i, n +i 
be uniquely determined irrespective of the order of increase of the indices m and n. 
That is, 

UVS„V~ l U~ l = VUS^U-W- 1 (13) 


which implies 

[U- 1 V- 1 UV,S i ., m>n \= 0 (14) 

For this to hold for all generic S t we must require U^ 1 V~ 1 UV = uj to be proportional 
to the identity matrix. Clearly uj satisfies det(u;) = u n = 1, so we obtain 


UV = ujVU , u = e i2 < 


(15) 


with v an integer 0 < v < n. U and V then satisfy Weyl’s braiding condition which 
characterizes a noncommutative (‘quantum’) torus |30| . 

The spin matrices Si ;m>n are now expressed as 


Si- m , n = u m v n s i v- n u- m = v n u m s l u~ m v- 


(16) 


Inserting the expressions (pT| , pT|) in the hamiltonian (]T|) we obtain the hamiltonian 
of the reduced system. As usual, the resulting hamiltonian has an infinite factor, 
due to the summation of the hamiltonians of the infinitely many identical cells over 
the complex plane. Dropping this trivial infinity, the reduced hamiltonian includes 
the kinetic terms of the fundamental cell and the interaction potential of particles 
in this cell with all other particles in all cells: 


^1 2 1 v ^ tr (U m V n S i y- n U- m S j ) 

- 2 Pi + 2 TO) ^1 00 (Xij + mci + nc 2 y 

where we adopted the notation x^ = Xi — x r 


( 17 ) 
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3 Noncommutative spin interaction potentials 

To proceed, we must identify the possible forms of U , V. We need the irreducible 
representations of the relation Call k the greatest common divisor of v and n. 
Then n = km and v — kfj,, for relatively prime m, //. The irreducible representations 
for U, V are m-dimensional ‘clock’ and ‘shift’ matrices. By a global U( n) spin trans¬ 
formation we can diagonalize either of U, V. Choosing U diagonal, the general form 
of U and V will be the direct sum of k of the above irreducible representations: 

U — diag{e i(t>0 , ... e^ fc_1 } ® u , V = diag{e i6 °,... e^- 1 } ® v (18) 

where (f> q , 9 q are arbitrary phases, determining the Casimirs U m and V m , and u, v 
are the m-dimensional clock and shift matrices 


u a p = S a/ s , v a j 3 = 5 a +i,p (mod m) , a, (3 = 0,... nr — 1 (19) 


So the acceptable U and V depend on 2k arbitrary parameters. 

To take advantage of the form (|T8|) for U, V we partition Si into k 2 blocks of 
dimension m x m each by using the double index notation 


{Si)ab = (SiYJp , a = pm + a , b = qm + /3 
The U(n) Poisson brackets in this notation are 

{(Si)3, &)%} = -ld l3 [(5i)S <57/3 <5rg d rq (S)# 
The m, n-surns that appear in (|T7|) then become 


& 


\QP 


,—im(t>p q —in0p q 


J S,) ° +n ^ n WJ ' A “ (*« + m °! + 


( 20 ) 


( 21 ) 


( 22 ) 


where the term rn — n — 0 is omitted if i — j. 

The above gives a potential interaction between particles i and j in the form 
of a modular function in Xij which depends on the spin components of particles i 
and j. To make the noncommutative character of the spin interaction manifest, we 
perform a change of basis in the spin states and define 

(23) 

a 

This is essentially a discrete Fourier transform in the sum of the a, [3 indices of S^. 
(Note that, for m odd, is actually antiperiodic in the index a if (3 is odd, and 
vice versa. Although we could have defined a properly periodic matrix, we prefer 
this slight inconvenience in order to make the ensuing formulae more symmetric.) 
In fact, it will be convenient to assemble the double indices (ex, (3) and (■ m,n ) into 
vectors ex and m. Similarly, we define c = (ci, c 2 ) and = (0 p , 6 p ). 
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The Poisson brackets of the Si are found from Q2T]) 


ax(3 


{(Si)?, (SjYfl = ^ S rs (s^ +f - ui ^ s r , (s^v 


(24) 


This is a structure extending the Moyal (star-commutator) algebra, the exponent of 
u) being the cross product of the discrete ‘momenta’ a and (3. For (rs) = (pq), in 
particular, it becomes the torus Fourier transform of the Moyal bracket 


{(Si)?, &)%} = HiM - ^) b X Pi (Si)? +(J 


(25) 


where 


0)2 


U) 2 


[x\ u i _j_ (26) 

0)2—0) 2 

is the o)-deformation of x. This is the so-called trigonometric algebra with periodic 
discrete indices m- 

Finally, by inverting (^3|) and substituting in (|22|) , the potential energy W in 
terms of the Si acquires the form 


"'-EE(S)? (§,)% 

i,j a-,p,q 


(27) 


The above includes two-body interactions, for i ^ j, as well as spin self-couplings, 
for i = j, arising from the interaction of each particle with its own images in different 
cells. The two-body potential W~ q (x) is 


1 _ , ixm p i(j> P q-m 

WZ^(x) = — N -— 

m _ (x + c • m) 2 

m \ / 


(28) 


while the spin self-coupling WZ q is 


1 . ,axm e i(j)pq-rn 

wz q = — y - _-_ 

m (c • m ) 2 


(29) 


If the above potentials were independent of the U(n) indices a and p, q, the sum 
over U{ n) indices in the potential energy expression ([27]) would simply be a U(n) 
trace and would give the U (n)-invariant coupling between the spins of particles i 
and j multiplying the standard Weierstrass potential of the elliptic Calogero model. 
In the present case, however, the above potential is spin-dependent and breaks U(n) 
invariance, introducing a star-product twist in the indices a and phase shifts (j) p in 
the indices p, q. Generically, the U{ n) invariance of the original model is broken 
down to an abelian U( l) fc , amounting to the transformation 


(S^ 


a/3 


o %L Pp 


(Si) 


\P1 p-Wq 


af3 


(30) 
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If (f) p are equal for k' values of p, the remaining symmetry U(l) k ' is enhanced to 
U(k'), corresponding to mixing the corresponding /^-components. 

The case u = 1, <j) p = 0 reduces to the standard spin-elliptic Calogero-Moser 
model. The case m = 1 (and thus u = 1) reproduces the U (n)-noninvariant model 
introduced in . The general case with uo ^ 1 is a new classical integrable model of 
the spin-Calogero type with a spin-dependent potential which is a modular function 
of the two-body distance ay,-. 

The sums appearing in 


and (|29|) could in principle have ambiguities due 
to the logarithmic divergence of the summation over the radial coordinate on the 
complex plane. This is, indeed, the case for the standard Weierstrass function and 
a specific prescription is needed to regularize it. Different prescriptions lead to 
different additive constants in the result. In our case, however, the presence of the 
extra phases renders the sums convergent and there is no regularization ambiguity. 

The potentials can be expressed in terms of theta-functions. WZ g (x) is a modular 
function on the complex torus (ci,c 2 ) with quasiperiodicity 


Wj?(x + Cl ) = e 
Wf 9 (x + c 2 ) = e 


_ c-typq + ^ra 2 


—ic> — 2m 


Oil 


wTM 

WP(x) 


(31) 


It has a double pole at x = 0, with principal part 




i 

mi 2 


+ 0 ( 1 °) 


(32) 


and no other poles in each cell. These properties uniquely define W$ q (x) and allow 
for an expression in terms of theta-functions. We put 


WP(x) = A uVt 


$1 1 

feU-O.f 

) $1 1 


d, ( 


2 


(33) 


where Q \ >2 are the as yet unknown zeros of W^ q {x) and the theta-functions appearing 
above have complex period T = c 2 /ci. This has the right quasiperiodicity under 
x —> x + Ci. In order to also have the right quasiperiodicity under x —> x + c 2 , Q i >2 
must satisfy 

Qi + V 2 — — x - 1 - ol ■ c (34) 

2tt m 

and to have the right behavior around x = 0 we must further have 


^ (jgO , ^ feg4 

* (so.) * 


• 0a6 

l - 

7T 


2i—a\ 

m 


A 


7T 2 ^(o) 2 

mc?M?)Mf) 


(35) 

(36) 
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The equations ([FI|) and (|55| ) above determine Q i and Q 2 , while ([SB]) in turn deter¬ 
mines A. The self-coupling W pq can then be extracted from W$ q (x) as 


W p d q = lim W p A(x) - 


x —>0 


mar 


(37) 


The sums appearing in (|28D and (|29|) are in general convergent, due to the pres¬ 
ence of the phases. For uj — 1, however, the phases are absent and terms with p = q 
have an additive ambiguity due to the need for regularization for the expression 
(^). In the theta-function expression this manifests in the fact that the equations 
for Q i i2 (|34| , |35D are satisfied for any Q i = — Q 2 . By applying the addition formula 


$ 1(2 + Q)til( x - <?)$4(0) 2 = $l(U 2 1?4(0 2 - J?4(^) 2 ^l(C) 2 


(38) 


this is seen indeed to amount to an arbitrary additive constant to the expression 
for W pp (x). The same holds for terms p, q for which rf> pq = 0. Such arbitrariness, 
however, corresponds to trivial redefinitions of the model by addition of constants 


of motion. This is explained in [22] and will not be elaborated here 


4 Conclusions and open questions 


In conclusion, we identified an integrable generalization of the elliptic spin model 
which breaks the spin U(n) invariance and promotes the potential to a modular func¬ 
tion introducing noncommutative spin twists. (Quantum generalizations of elliptic 
spin models have appeared recently, but with U( n) invariant interactions p2| .) 

There are clearly many issues that deserve further study. The conserved quanti¬ 
ties and Lax matrix of this model can, in principle, be obtained from the correspond¬ 


ing quantities of the unreduced model; this was done, e.g., in [gg] for a specific case. 
A derivation in the present case would be useful. Further, the modular potential 
involves implicitly defined Q 1 ^; a more explicit and symmetric expression would be 
desirable. 

The properties of the spin interaction should also be clarified. I 11 particular, 
it would be interesting to see if some deformation of U( m) can be identified as a 
dynamical symmetry. 

Finally, the quantum mechanical extension of the model is, perhaps, the most 
interesting and pressing question. This is not trivially accessible by the method 
used here since, in general, the constraints implied by the phase space restrictions 
are second class and we cannot carry over the solution of the unrestricted quantum 
model and apply the constraints as operator relations on the Hilbert space. The 
above issues remain interesting topics for future investigation. 

Acknowledgement: I would like to thank C. Zachos for useful comments on the 
manuscript. 
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